Topological signatures of globular polymers 
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Simulations in which a globular ring polymer with delocalized knots is separated in two interacting 
loops by a slipping link, or in two non-interacting globuli by a wall with a hole, show how the minimal 
crossing number of the knots controls the equilibrium statistics. With slipping link the ring length 
is divided between the loops according to a simple law, but with unexpectedly large fluctuations. 
These are suppressed only for unknotted loops, whose length distribution shows always a fast power 
law decay. We also discover and explain a topological effect interfering with that of surface tension 
in the globule translocation through a membrane nanopore. 

PACS numbers: 36.20.Ey,02.10.Kn, 87.15.Aa 
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Most often, ring polymers are experimentally studied 
in situations in which their topological entanglement does 
not change in time. However, most of our theoretical and 
numerical understanding of polymer statistics relies on 
ensemble descriptions in which the rings assume all possi- 
ble topologies 'ill . An open challenge is that of determin- 
ing up to what extent specific permanent entanglements 
in the form of knots or links Q can affect thermodynamic 
quantities and what is their possible role in determining 
peculiar behaviors when the polymer is subject to geo- 
metrical constraints interfering with the topology. 

Different topologies are expected to determine differ- 
ent non-extensive corrections to the free energy of a sin- 
gle ring in the limit of infinitely long chain. For exam- 
ple, the prime knot components of ring polymers in good 
solvent are weakly localized in this limit [Sl, |j] . As a con- 
sequence, if N is the chain length, each component de- 
termines a correction ~ fc^ hi(N)/N to the entropy per 
monomer 0, @| . Indeed, each component behaves asymp- 
totically as a point-like decoration which can place itself 
anywhere along the ring. Radically different conditions 
are realized below the theta temperature Tg Q. Indeed, 
in the globular phase knots are expected to delocalize. 
Numerical simulations indicate that the topological en- 
tanglement spreads on average on a portion of the ring 
whose length is proportional to N [1, 0]. Fixed topol- 
ogy is also known to determine a finite size correction 
to the free energy of a globular ring, which is asymptot- 
ically negligible in comparison with that due to surface 
tension (~ N~-^/^). Indeed, at a temperature T ^ |Tfl, 
the correction per monomer has been estimated [10l | as 
^ Cn'^N~'^ where a ~ 1.45, ric is the minimal number 
of crossings of the knot 0| , and C is a remarkably large 
negative amplitude. Thus, ric qualifies as a topological 
invariant possibly relevant for the thermodynamics of a 
globular ring. 

In the present work we face the challenge of elucidating 
the role of this invariant, by establishing some empirical 
laws through which it rules the statistics of the globule 
when suitable local geometrical constraints are imposed. 



These laws, which generally hold in the presence of re- 
markably large fluctuations, give a precise meaning to 
the notion of knot delocalization. We also discover novel 
thermodynamic phenomenona which are a peculiar con- 
sequence of topology. 

Suppose one forces a ring polymer to pass into a slip- 
ping link which divides it into two loops. The link is 
narrow enough to prevent the passage of topological en- 
tanglement from one loop to the other. For instance, 
if the ring has a composite knot with two prime trefoil 
(3i) components 0, things can be arranged so that each 
fluctuating loop encloses one trefoil knot. The loops are 
kept unlinked. We model flexible ring configurations as 
TV-step self avoiding polygons on cubic lattice [ll|. An 
attractive interaction J between nearest neighbor visited 
sites allows to obtain a collapsed globular phase for low 
enough temperature T 0] (Fig. HI . A Monte Carlo simu- 
lation method adequate to preserve the ring topology is 
the grand-canonical BFACF one ill)]. However, staight- 
forward application of BFACF to the globule meets a dif- 
ficulty. Indeed, indicating by K the step fugacity and by 
Ztv the canonical partition function, the grand canonical 
average (N) = Y^nNK^ Zn /^nK^ Zn does not grow 
continuously to +oo upon approaching from below the 
critical value Kc of the fugacity. To the contrary, one 
gets evidence of a discontinuous infinite jump right at 
K = Kc [1]. We understand here this behavior in the 
light of the expected ■ [l^ asymptotic form of Zj^: 



MZn) = -FN/kT - const + fiN + aN^'^ -f 
(a-2)ln(A^) + ^ 



(1) 



where ji = — In(ifc), cr < is the interfacial tension, 
and a is an unknown specific heat exponent that is ex- 
pected to be independent of topology. The last term on 
the r.h.s. of Eq.(I2|) is the above mentioned topological 
correction to the total free energy Fm- The presence of 
the surface term aN^^^ implies that (N) can not di- 
verge continuously to +oo for K approaching Kc from 
below. In order to allow a continuous growth of {N), 



FIG. 1: A'^ = 2040 globule with slip link separating a 3i knot 
in loop 1 (yellow) from a 4i knot in loop 2 (blue). 



we choose to multiply the usual grand-canonical weight 
of the ring configurations by a factor exp[(Af — N^)'^ /v] 
which forces TV itself to fluctuate around a value close to 
A'o if V is chosen small enough. Using still BFACF moves 
the simulation augmented with the new statistical weight 
allows a quasi-canonical sampling of configurations with 
N close to the value determined also by Nq and v. To 
improve sampling efficiency, we also implement a multi- 
ple Markov chain scheme 13| over different values of A^o 
keeping K fixed. Nevertheless, the simulations need a 
long CPU time (months) to sample a consistent statis- 
tics. For this reason we choose to sample only at the 
temperature T — 2.5 J, as in [lo| . 

Extensive simulations of the 3i vs 3i configuration al- 
low to sample for various restricted ranges of N the prob- 
ability density function (PDF) of h/N, P(/i/iV), where 
li is the fluctuating length of one of the loops (Fig. [2^ 
). Remarkably, even for large N the histogram does not 
seem to present the bimodal shape found in the good 
solvent case and indicating a dominance of configurations 
with large unbalance between li and I2 [Mill- To the con- 
trary, an asymptotically flat histogram seems compatible 
with the data. Thus, while the two loops on average 
share equally the total ring length (Fig. 2b), there are 
very broad fluctuations because all possible partitions of 
the total ring length are almost equally probable. 

To gain further insight, we simulate also the case in 
which the loops are both unknotted. In this case strongly 
unbalanced situations are clearly fovoured. The PDF of 
the length of the smaller loop, say h, shows a power 
law decay - l^'', with x = 1.55 ± 0.04 (Fig. E]). This 
suggests that the metric exponent of this loop is = 
a;/d«0.52. A value of ly close to 1/2 is consistent with 
the expectation of Gaussianity of a collapsed chain on 
relatively small length scales [IJ] . When for example only 
loop 2 has a 3i knot, while loop 1 is unknotted (0), we 
observe that li never grows substantially compared to iV, 
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FIG. 2: (a) Histograms of P{li/N) for 3i vs 3i. Different 
curves correspond to different A'^ values (see legend), (b) 
{Ii)n (circles) and {12)n (squares) as a function of A'', (c) His- 
tograms of P{li/N) for 3i vs 4i and (d) corresponding aver- 
age loop lengths, (e) Histograms of P{li/N) for 6i (loop 1) 
vs 3i#3i and (f) corresponding average loop lengths. 



and the same power law behavior holds for the PDF of li 
(Fig. O. This power-law, implying that (h) - N°-^^^^\ 
denotes a weak localization of the unknotted loop. 

The above results are fully consistent with delocaliza- 
tion of the 3i knot inside the ring, since in no case a loop 
containing the knot displays a stable regime in which its 
average length is a vanishing fraction of N. Similar re- 
sults hold if 3i is replaced by another prime knot. We 
further examine a competition 3i vs 4i. The plots re- 
ported in Fig. [21; show that also in this case the lengths 
of the loops keep fluctuating very broadly for increas- 
ing N, while P(li/N) is not symmetric with respect to 
h/N = 1/2 anymore. Quite remarkably, in this and sim- 
ilar competitions (3i vs 7i, 4i vs 61 etc.), a simple law is 
well obeyed by the canonical averages {Ii)n and {12)n- 
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FIG. 3: Log- log plots of P{li) for vs (upper curves, for 
three values of A*'), and vs 3i (lower curves, shifted one 
decade down). The dashed line represents a power-law l^^'^. 
Each A'^ includes the statistics in the interval [A'' — 50, A'' -|- 
50]: this does not alter the power-law tail and increases the 
statistics considerably. 




FIG. 4: Example of a wall-globule system. The hole (white) 
separates the A'' = 1000 chain into two globuli with 3i knot 
(colored differently). 



In this case the loop with the 3i knot obtains on av- 
erage a fraction of the chain length equal to 3/7 w 0.43, 
while 4/7 « 0.56 go to the loop with 4i knot. This is 
evidenced in Fig. by the linear fits of the average loop 
lengths as a function of N. 

The law in Eq.(2) highlights the key role played by the 
topological invariant ric in the delocalized regime. This 
role is further enphasized by considering the competition 
of two knots that are different but with the same Uc. 
Fig. [2^ and [2f show the results for 3i#3i vs 6i. We see 
that, while the fluctuations remain very broad, also in 
this case the (/i)jv's still obey Eq. In addition the 

shape of the histograms is quite symmetric as in the 3i vs 
3i case. However, the possible convergence towards a flat 
histogram is definitely slower than in that case. These 
results suggest that the number of prime components of 
a knot is not a relevant invariant, unlike in the swollen 
regime. 

In experiments the slipping link could be, e.g., a short 
(unknotted) ring which is not linked to the knotted one. 
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FIG. 5: Histograms of P{li/N) for 3i#3i (loop 1) vs 6i (top 
panel) and for 3i (loop 1) vs 4i (bottom panel) when the 
two globular loops are separated by an impenetrable wall. 
Different curves correspond to different N values. 



but just constrains it to pass through its interior. To 
discover further consequences the topology may have in 
the globular phase we replace the slipping link by a suffi- 
ciently narrow hole in an impenetrable wall. Fig. |4l This 
schematizes the translocation of a ring polymer through 
a membrane or solid state nanopore, a phenomenon that, 
thanks to recent progress in nanotechnology can be in- 
vestigated experimentally 15, [l^. A related biological 
issue is the effect that knots may have on the election 
process of packed viral DNA into the host cell 17, 18| . 
Due to the presence of the wall the two loops do not 
interact anymore and can be considered as two indepen- 
dent, knotted globuli in competition. As shown in Fig. [S] 
for a case with two different knots having however the 
same rid = = 6, P{li/N) develops two remarkably 
symmetric peaks for N > 1000, which become separated 
by a very high free energy barrier already at N ^ 2300, 
the maximum N for which our multiple chain sampling 
works efficiently. Similar results are obtained for the sim- 
pler case 3i vs 3i. 

The situation turns out to be totally different when 
"■ci 7^ nc2- In these cases we observe an asymmetry, 
indicating a clear dominance of the globular loop hosting 
the knot with higher Uc- This dominance becomes more 
and more pronounced as N increases. We report the 
results for the case 3i vs 4i with rings of length up to 
N = 4000 in FigEl 

How topology can produce the asymmetry manifested 
by the above results for the wall case can be explained as 
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follows. The equilibrium share of the ring length between 
the two globuli should stably minimize the free energy. 
Due to independence, the total free energy is simply the 
sum of the free energies of the two globuli. We can as- 
sume that in the relevant ranges explored by our simu- 
lations both li and I2 are sufficiently large for Eq.([2]) to 
give a reasonable approximation of the globule free en- 
ergy. The surface term in Fn, if taken into account as 
the only correction, implies that both configurations in 
which one of the globuli has length < N, while the other 
is close to minimal, realize equivalent stable minima of 
the total free energy with constraint h + h — N. Let 
us indicate by mi and m2 the value of h and l2, respec- 
tively, in the stable configuration in which it happens to 
be the shorter loop. From Eq.Q one gets for the ratio 
of canonical probabilities of these configurations: 



In 



P{h -1712,12 ^ma) 
P{li = TOi,/2 -mi) 



{N - m2)2/3 + ml^' -{N- mif'^ - 



2/3 



{a - 2) In 



{N - m2)m2 
{N — toi)toi 



C 



N — 7712 



m2 



mi N — mi 



(3) 



In the case rid — 't-c2 the presence of the topological cor- 
rection does not matter, and thus we expect to see with 
the same frequency a dominance of either one of the two 
globuli. The only topological effect left could be due to a 
slight difference between mi and m2 in the cases in which 
the two knots in 1 and 2 are not identical. However, this 
difference should be extremely small. Indeed, for exam- 
ple the minimal lengths of lattice knots with the same ric, 
which are lower bounds for mi and m2 , are known to be 
almost coincident [l9| . When we instead have rid 7^ nc2, 
the two configurations favored by surface tension are not 
equally probable anymore. The term apt to determine a 
substantial deviation of the quantity in Eq.(3) from zero 
originates in fact from the topological correction: 



C 



"■c2 
TO2 



(4) 



The main reason for this deviation is the remarkably large 
absolute value of C, which amplifies any even slight dif- 
ference — To get an estimate of the quantity in 
Eq.dH) for the first knots, we calculated it for C « —170 
and a ~ 1.45, the values found for T/ J = 2.5 [l3|, and 
mi and m2 replaced by the existing estimates of the min- 
imal lattice knot lengths [l^ , consistently with the peak 
structures reported in FiglS] For the case 3i vs 4i, the 
value of (jU turns out to be ~ 7.5. 

Thus, a topological mechanism explains why a stable 
configuration is that in which the globule with smaller 
ric is reduced to a minimal size, while the other takes 



most of the ring length. This novel phenomenon should 
be relevant, e.g., for a slow dynamics of translocation of 
globular knotted ring polymers through membrane pores. 

Putting things in perspective, althought necessarily 
limited to the less complex knots, our results are accurate 
enough to clearly qualify ric as the invariant controlling 
delocalization in the globular phase and determining un- 
expected thermodynamic effects of genuinely topological 
nature. The outlined scenario is totally different from 
that expected in the swollen polymer regime. Even if un- 
usually large chain lengths need to be explored in order to 
fully clarify their asymptotics, the discovered phenomena 
are most relevant in practice for finite N . 
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